function diff = int_delta_fun(m,param)

    diff = nan(size(m));
    f = @(x) ( (param.zeta+delta_fun(x,param))./x);

    for i=1:length(m)  
        diff(i)=integral(f,param.m_l,m(i));
    end

end